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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

DRAG DUE TO LIFT AT MACH NUMBERS UP TO 2 .0 
By Edward C . Polhamus 


INTRODUCTION 


In references 1 and 2 it has been shown that, if the "area rule" is 
utilized properly, it is possible to obtain values of zero-lift drag which 
for a wide variety of wing- fuselage configurations, approach that for the 
basic fuselage alone. This fact malces the selection of a wing less depend 
ent on its zero-lift drag and therefore allows a wider range of wings to 
be considered with regard to drag due to lift. The purpose of this paper 
therefore is to discuss the effect of wing geometry on the drag due to 
lift at Mach nvunbers up to 2.0 and the effect, of application of the area 
rule on the drag at lifting conditions . 

In figure 1 a typical variation of the drag with lift coefficient 
for a plane, or flat, wing is shown by the solid line on the left-hand 
side of the figure. For a plane wing the minimum drag occurs at zero 
lift and theoretically has a parabolic shape with the increment due to 
lift equal to a constant times the lift coefficient squared. In 

general, the data for the wings presented in this paper were fairly 
linear plotted against up to lift coefficients of about 0.5 and 

therefore the slope will be used to describe the drag-due-to- 

lift characteristics of plane wings in this lift range. For a cambered 
or cambered and twisted wing the drag curve, as shown by the dashed line, 
does not have its minimum at zero lift and therefore the drag polars will 
be used to describe the chairacteristics of this type of wing. 

Now, if viscous forces are neglected, the drag due to the lift can 
be divided into two components - a thriist conponent of the suction force 
caused by the flow about the nose of the airfoil, and a drag conponent 
of the normal force. For a two-dimensional wing these two components 
exactly balance each other; however, for a three-dimensional wing the 
drag conponent of the normal force is greater than the thrust component 
of the suction force, since a higher angle of attack is required to 
develop the same lift, and an induced drag results. At subsonic speeds 
the rate of change of the induced drag with lift squared can be approxi- 
mated by l/jtA as illustrated in the bottom part of the right-hand side 
of figure 1. Additional drag also occvirs if the suction force is not 
fully developed at the leading edge. For the extreme case of zero suc- 
tion the drag due to lift is equal to the conponent of the normal force, 
and the rate of change is therefore equal to the reciprocal of the 
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lift-curve slope as illustrated in the top part of the figure. The 
drag curve of a wing usually lies somewhere between these two extremes 
and its relative position between these two limits is dependent to a 
large extent on the amount of suction developed at the leading edge and 
is therefore a function of such parameters as Reynolds number, Mach num- 
ber, thickness, and leading-edge radius. The two limits, of course, are 
primarily a function of plan form and Mach number. 


EFFECT OF REYNOLDS NUMBER 


Figure 2 shows the effect of Reynolds number on the drag due to 
lift of an aspect-ratio-2 delta wing having an NACA OOO 5-65 airfoil 
section (ref. 5 and unpublished data). The results are presented in 

the form of the drag-rise parameter against Reynolds number 

for several Mach nimbers. Also shown are the subsonic and the M = I .7 
theories for full leading-edge suction and the values for zero suction 
given by 

The results indicate that at a Mach number of 0.25 there is a 
rather large increase in drag due to lift with decreasing Reynolds num- 
ber but that as the Mach number increases the effect of Reynolds number 
diminishes and is relatively unimportant at a Mach number of 1 . 7 . The 
increase with decreasing Reynolds nvuriber is probably due in part to the 
fact that the combination of low Reynolds number and a relatively sharp 
leading edge is conducive to leading-edge separation resulting in a loss 
of leading-edge suction. In addition, a part of this variation is 
probably due to the fact that at low Reynolds numbers the transition 
point moves forward with increasing lift resulting in an increase in 
viscous forces with lift. The decreasing effect of Reynolds number 
with increasing Mach numbe.r is due to the fact that the difference 
between the theory and the zero-suction case decreases . with increasing 
Mach number and the fact that the flow about the leading edge is 
affected by compressibility. It should be pointed out that, while the 
Reynolds number based on the mean aerodynamic chord was used here to 
define more clearly the variation with Reynolds number for a given wing 
it appears that the drag due to lift at a given Mach nuiriber is more * 
dependent upon the Reynolds number based on leading-edge radius. A 
recent correlation (ref. 14-) based on this parameter succeeded in bringing 
the drag-due-to-lift parameter into fair agreement for a large number of 
aspect-ratio-2 delta wings having various airfoil sections. It should 
also be pointed out, however, that, for plan forms where compressibility 
®ff®cts are a function of thickness ratio or leading— edge radius, corre- 
lations based on the leading-edge Reynolds number would not be expected 
to bring the data into agreement at all Mach numbers . 
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EFFECT OF THICKNESS 


Figure 5 illustrates the effect of wing thickness ratio on the drag- 
due-to-lift factor for unswept wings of aspect ratio 4 at Reynolds num- 
bers of approximately U x 10^ (refs. 6 , and unpublished data from the 
Langley l6-foot transonic tunnel) . In addition to the experimental data, 
the theory for full suction is also shown. It will be noted that at sub- 
sonic speeds a decrease in thickness ratio from 8 percent to 4 percent 
increased the drag-due-to-lift factor; for example, at a Mach number 
of 0.6 it was increased by approximately 6o percent. This increase with 
decreasing thickness ratio is probably due to the fact that the h-percent- 
thick airfoil section has a considerably smaller leading-edge radius and 
therefore develops less leading-edge suction. However, it will be noted 
that as the Mach number is increased the c\irves tend to converge and at 
a Mach nvuriber of about 0.88 there is little effect of thickness. This 
is due to the fact that, although the thick wing develops more suction 
at low speeds, the effect of compressibility on the flow about the 
leading edge is greater than for the thin wing. Above a Mach nvimber 
of 0.88, the 4-percent- thick wing has considerably less drag due to lift 
than the 6-percent- and 8-percent -thick wings due to the fact that in 
this Mach number range the resultant force is normal to the wing chord, 
and since the thin wing has the higher lift-ciorve slope it has the lower 
drag due to lift. This is illustrated by the two dashed curves repre- 
senting the reciprocal of the lift-curve slope for the 4-percent- 
6-percent-thick wings. 

Figure 4 shows the effect of thickness on the drag due to lift of 

a delta wing of aspect ratio 2 at a Reynolds nxmiber of 5 x 10^ (ref. 5) . 

At subsonic Mach numbers it will be noted that the results are similar 
to those for the unswept wings (fig. 5) with the thin wing having the 
highest value. of drag due to lift. However, at the higher Mach numbers 
the effect of thickness did not reverse for the delta wing’ as it did 
for the unswept wing and the thin wing still had the highest drag due 
to lift. It will also be noted that even at the highest Mach mmiber 
tested the drag is lower than the reciprocal of the lift-curve slope, 
an indication of some suction being developed. This is due to the fact 
that the Mach number normal to the leading edge of this wing never 
exceeded a value of about 0.80. The vertical dashed line represents 
the free-stream Mach niimber for which the Mach number normal to the 
leading edge is equal to 0.9 which is approximately eq\ial to the Mach 
number of the unswept wings for the case of zero suction. In order to 
indicate the variation with Mach number in the transonic range, the 
results of a rocket-propelled model of similar plan form having a 

thickness of 6^ percent (ref. 7) is shown by the long and short dashed 
curve. 
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EFFECT OF LEADING -EDGE RADIUS 


The effect of leading-edge radius on the drag due to lift of an 
unswept wing (ref. 5) is illustrated in figure 5- The wing had an 
aspect ratio of 5 ^ a- taper ratio of 0 . 59 ^ and a thickness of 5 per- 
cent and was tested with a biconvex section and with a biconvex sec- 
tion modified with an elliptical nose having a radius of 0.045 percent 
of the chord. It will be noted that the results are similar to those 
obtained in the thickness investigation, with improvements with 
increasing leading-edge radius occurring only at subsonic speeds. It 
should be pointed out that the two curves are coincident at supersonic 
speeds . 

Figure 6 presents the results obtained on a 45° swept wing of 
aspect ratio 4 which was tested with several modifications to the basic 
NACA 65 AOO 6 airfoil section in the Langley high-speed 7- 4»y 10-foot 
tumel. The three configurations tested were a sharp edge having zero 
radius, the normal radius of 0.24 percent chord, and a radius of 
0.72 percent chord. Inasmuch as only a limited Mach number range was 
covered in these tests, the results are presented as plotted 

against Cj^ at a Mach nvimber of O. 9 O. The results indicate that no 

improvement occirrred with increase in the leading-edge radius at this 
Mach nxamber. 


EFFECT OF ASPECT RATIO 


Flgvire 7 illustrates the effect of aspect ratio on the drag due to 
lift through the Mach number range. The wings were of delta plan form 
and 5 percent thick and had aspect ratios of 2 and 4 (ref. 5 ) • The 
results indicate, as would be expected, that the higher aspect ratio 
has the lower drag due to lift throughout the Mach number range. How- 
ever, it will be noted that the difference between the two aspect ratios 
is considerably greater than that indicated by the theory. It will be 
noted, however, that the effect of aspect ratio on the reciprocal of the 
lift-curve slope, which represents the zero-suction case, is approximately 
twice that for the full-suction theory at subsonic speeds . The larger 
effect of aspect ratio obtained in the experiments is therefore not sur- 
prising since these thin wings lose a good portion of suction. 


EFFECT OF SURFACE SHAPE 


The previous figures have illustrated the effect of various parameters 
on the drag due to lift of planar wings and have shown that, in general. 
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the drag due to lift is considerably higher than the theoretical values 
largely because of separation at the nose and the accompanying loss of 
thrust. However, a theoretical study by Jones (ref. 8) has shown that 
an effective leading-edge thrust can be obtained by cambering and 
twisting the wing. 

Figure 8 presents the results of cambering and twisting a swept 
wing of aspect ratio 4 (ref. 9 and vinpublished data from the Langley 
8-foot transonic tunnel and the low- turbulence pressure tunnel) . On the 
left-hand side of the figure the results are presented for the case of 
the cambered and twisted wing having incidence at the fuselage 

which results in low fuselage angles at moderate lift coefficients. The 
results are presented as a plot of Cp against Cj^ s-t a Mach number 

of 0.9 for the plane wing at zero incidence, for the wing cambered and 
twisted for a uniform load distribution at a Cl of O.il- and M = 1.2 

(e = 4.5°), and for the wing cambered and twisted for a triangixlar span 
load and a rectangular chord load at a Cl of 6'’.4 and M = 0.9 (e = 15°) . 

The triangular span load of the latter case was used in an attempt to 
in5)rove the pitching -moment characteristics. The results indicate large 
increases in drag for both the camber and twist distributions . However, 
on the right-hand side of the figure, results are presented for the wing 
cambered and twisted for a uniform load tested on a slightly different 
fuselage but having approximately zero incidence at the fuselage. These 
results indicate substantial reductions in drag above a lift coefficient 
of about 0.15 for the cambered and twisted wing. For the case of zero 
Incidence the fuselage is developing lift at the design condition and 
therefore the wing-fuselage combination represents the wing alone for 
which the camber and twist were designed considerably better than the con- 
figuration having 4.5° incidence which results in low fuselage angles 
in the moderate lift range. 


Figure 9 presents the^results of an aspect-ratlo-2 delta wing tested 
with three different surface shapes (ref. 5) ‘ a planar siarface, a sur- 
face cambered and twisted for a trapezoidal spanwise load distribution, 
and a surface which was planar over the inboard 80 percent of the local 
semispans. This third siarface was a modification of the surface required 
for an elliptical span loading and was used in order to sin^jlify con- 
struction. It should be mentioned that the wing incidence was zero at 
the plane of symmetry for all three cases. At a Mach number of O.91, 
it will be noted that both types of surface modification resulted in 
improvements in the drag characteristics but that the sinple nose camber 
was superior below a lift coefficient of about 0.5. As the Mach number 
increased the improvements diminished for both surfaces and at a Mach 
number of I.55 no improvement was obtained; however, less penalty occurred 
for the simple nose camber. 
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Figvire 10 shows the results obtained from tests in the Langley 
8-foot transonic tunnel of a similar wing in which the extent of the 
nose camber was varied. The wing had an aspect ratio of 2.2 and a 
modified MCA 0004-65 airfoil section. Two nose cambers were tested 
and both were of constant chord, one being 4 percent (modification A) 
and the other 8 percent (modification B) of the mean aerodynamic chord. 
The camber covering 4 percent of the mean aerodynamic chord was obtained 
by shearing the ordinates so that the bottom surface was parallel to the 
chord line . The camber covering 8 -percent of the mean aerodynamic chord 
was obtained by extending the chord and displacing the leading edge an 
amount equal to 1.3 percent of the longitudinal distance X from the 
wing apex. This was a modification of the surface shape required for 
an elliptical loading at a lift coefficient of 0.15 and was similar to 
that presented in figure 9 except that it was of constant chord. On 
the left-hand side of the figure drag polars are shown for the basic 
wing and the two modifications at a Mach number of 1.0. The results 
indicate that modification A (4 percent) had no effect on the drag 
while modification B (8 percent) resulted in a substantial reduction in 
drag except for an extremely small increase at zero lift. Although the 
coefficients are based on the actual areas, it should be pointed out 
that even the actual drag for a given lift is less for modification B 
than for the basic wing. The effect of Mach number on the variation of 
drag at a lift coefficient of 0.3 for the three configurations is shown 
on the right-hand side of the figure. It will be noted that both modi- 
fications resulted in improvements at the lower Mach number but that 
modification A had no effect above a Mach number of about 0.90; however, 
modification B resulted in improvements throughout the Mach number 
range investigated. 


EFFECT OF TRIMMING 


In order to reduce the weight and zero-lift drag of an aircraft, 
tailless configurations are sometimes used. However, since a tailless 
design, in general, obtains its trim from a surface on the wing, large 
deflections of this surface are required because of the short moment 
arm. These large deflections, of course, result in additional drag 
which could have an important effect on the performance. This is 
especially true at supersonic speeds because of the increased stability 
caused by the rearward movement of the wing aerodynamic center in going 
from subsonic to supersonic speeds and the higher drag due to flap 
deflection. Figure 11 illustrates this effect of trimming on the varia- 
tion of the drag with lift. The model was an aspect-ratio-2 delta wing 
having an NACA OOO5-63 airfoil and a constant-chord flap equal to 10 per- 
cent of the wing mean aerodynamic chord (ref. 10) . At a Mach number 
of 0.90, it will be noted that a positive flap deflection of 4° resulted 
in a reduction in the drag due to lift. However, for a stable tailless 
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configuration negative deflections are required to trim the airplane 
through the positive lift range which results in an increase in the 
drag due to lift. At a Mach number of 1.90, the increase due to 
trimming the airplane is considerably greater than at O.9O because 
of the aforementioned increase in stability and drag due to flap 
deflection at supersonic speeds. 


EFFECT OF APPLICATION OF THE AREA RULE 


In references 1 and 2 it was shown that indentations of the fuse- 
lage according to the Mach number of 1.0 area rule resulted in large 
decreases in the zero-lift drag of wing-fuselage combinations at tran- 
sonic speeds. The question now arises as to whether these benefits are 
maintained under lifting conditions. Figure 12 shows the effect of a 
Mach number of 1.0 body indentation on the drag of a wing-fuselage com- 
bination. The wing had an aspect ratio of k, of sweep, a taper ratio 
of 0.5, and an NACA 65AOO6 airfoil section and was tested in the Langley 
8-foot transonic tunnel. The results are presented as total drag coeffi- 
cient against Mach number for lift coefficients of 0 and 0.5 and indicate 
that the large reductions in drag at transonic speeds due to body inden- 
tation were to a large extent maintained in the lifting condition. At 
supersonic speeds the Mach number 1.0 indentation had negligible effect 
at either lift coefficients of 0 or 0.5. 

Figure I3 shows the iiriprovement in the maximum lift-to-drag ratio 
associated with this application of the area-rule concept. The results 
for both the basic configuration and the configuration with the indented 
body are plotted against Mach nvunber and it will be noted that, below 
a Mach nvimber of about 1.1|-, the lift-to-drag ratios were improved and at 
a Mach number of 1.0 (the design condition) the increase amoionted to 
approximately 37 percent. 

At the present time little has been done in attempting to develop 
area distributions which might actually reduce the drag increment due 
to lift. However, figure l4 presents the results of one such investiga- 
tion conducted on a wing of aspect ratio 4 having 45° of sweep (ref. 11) . 
The basic body was cylindrical rearward of the wing leading edge and was 
modified by several types of indentations . The first indentation, 
designated by the letter @ in the figure, was symmetrical around the 
fuselage and was determined by the Mach number 1.0 area rule. The other 
two indentations tested were more abrupt indentations superiit5)osed first 
on the upper half (c) and then on the lower half @ of the symmetrical 
indentation. On the left-hand side of the figure the drag at zero lift 
is presented against Mach number and it will be noted that all the modi- 
fications gave about the same reduction in drag. On the right-hand part 
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of the figure the same comparison is made for a lift coefficient of 0 . 3 . 
The results indicate that the syTnmetrical area-rule indentation resulted 
in about the same reduction in drag as at zero lift which is consistent 
with figure 12. However, when the more abrupt indentations were added, 
additional reductions in drag resulted with the lowest occvurring for the 
indentation below the wing. 


COMBINED EFFECTS 


At transonic speeds it has been shown that application of the area 
rule and the use of camber and twist results in significant reductions 
in drag. Figure I 5 shows the effect of combining these two methods at 
transonic speeds. The tests were conducted in the Langley 8-foot tran- 
sonic tunnel on a model having sweep, an aspect ratio of a 

taper ratio of 0.6, and an NACA 65 AOO 6 airfoil section. The model was 
tested ( 1 ) with the basic wing and body, ( 2 ) with the basic wing and 
the body indented according to the area rule, and ( 3 ) with the wing 
cambered and twisted for a uniform load at C^ = 0.1<- and M = 1.2 in 

combination with the indented body. On the left-hand side of the fig- 
ure drag polars are presented for the three configurations at a Mach 
number of 1.0. It will be noted that indenting the fuselage resulted 
in large reductions in drag throughout the lift range. Camber and twist 
resulted in a rather large increase in minimum drag but resulted in 
improvements above a lift coefficient of about 0.2. On the right-hand 
side of the figure the maximum lift-to-drag ratios are plotted as a 
function of Mach number. At a Mach mmiber of 0.8 the iri 5 )rovement is due 
mainly to the camber and twist and resulted in an increase from I 3 to I 7 . 
At a Mach number of 1.0, the inprovement is due mainly to the body 
indentation and resulted in an increase from about 7.5 to 11 . 5 . 


CONCLUDING REMARKS 


In conclusion, it appears that Reynolds number has a rather large 
effect on the drag due to lift of thin wings at low speeds but that this 
effect decreases considerably with increasing Mach number. 

Comparisons of wings of various thicknesses indicate that at sub- 
sonic speeds an increase in thickness is beneficial, whereas, in general, 
at transonic and supersonic speeds no gains and possible losses occur 
imless the wing leading edge is highly swept which results in relatively 
low subsonic speeds normal to the leading edge. Similar results are 
indicated with regard to leading-edge radius . 
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Although camber and twist are effective in reducing the drag due 
to lift at the design condition, providing the correct wing incidence 
is used, it appears that simple nose camber will result in similar 
gains with less penalty near zero lift. 

The reductions in minimum drag associated with application of the 
area rule by means of fuselage indentations are maintained in the lifting 
condition and significant improvements in the lift-to-drag ratios result. 
In addition,' from preliminary tests, it appears that local modifications 
to the fuselage indentations may result in additional reductions in drag 
at lifting conditions. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 3^ 1953* 
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Figure 1 


EFFECT OF REYNOLDS NUMBER 
A*2i NACA 0005-63 
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Figure 2 
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EFFECT OF WING THICKNESS 
A=4;A«0*i X«0.5 



Figure 3 



Figure 4 
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EFFECT OF LEADING -EDGE RADIUS 
A = 3; A = 0°; X => 0.39; t/c = 0.03; R = 3.8X10® 



Figure 5 


EFFECT OF LEADING-EDGE RADIUS 
A=4jX*0.3; A*45<’j NACA 65A006 (BASIC) 
M=0.95 R=3 xI0® 
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Figure 6 
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EFFECT OF ASPECT RATIO 
t/c*0.03; R* 3 x 106 


ACq 



; Figure 7 

i 


EFFECT OF INCIDENCE ON CAMBER AND TWIST 
A = 4; A = 45": X = 0.6; NACA 65A006 






Figure 8 
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EFFECT OF SURFACE SHAPE 
A=2 ; NACA 0005-63 ; R = 7.5 x 10® 





Figure 9 


EFFECT OF EXTENT OF NOSE CAMBER 
A=2.2; NACA 0004-65 (MOD.)j R = 4x|o6 



Figure 10 
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EFFECT OF TRIM 


A=2: NACA 0005-63; R» 4.4X10® 



0 .1 .2 .3 .4 0 .1 .2 .3 .4 
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Figure 11 


APPLICATION OF AREA RULE 
A=4iA=45“; X=0.3; NACA65A006 



0 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 
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Figure 12 
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APPLICATION OF AREA RULE 

A=4;A=45°i X=0.3;NACA 65A006 



0 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 
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Figure 15 


EFFECT OF LOCAL BODY INDENTATIONS 

A=4;X= 0.6; NACA 65A006 
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Figure l4 
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APPLICATION OF AREA RULE AND CAMBER AND TWIST 



Figure 15 
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